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QCD sum rules are evaluated at finite nucleon densities and temperatures to 
determine the change of pole mass parameters for the lightest vector mesons p, 
uj and 4> in a strongly interacting medium at conditions relevant for the starting 
experiments at HADES. The role of the four-quark condensate is highlighted. 
A few estimates (within a fire ball model and BUU calculations) of dilepton 
spectra in heavy-ion collisions at 1 AGeV are presented. 

1 Introduction 

Hadrons are excitations of the QCD vacuum. There are many different species and 
masses of them. Their excitation spectra resemble, to some extent, to level schemes 
of atoms. The latter ones are known to be modified by external fields: The Stark 
and Zeeman effects manifest themselves as shifts of spectral lines. By analogy, one 
might expect a similar behavior of hadrons in strongly interacting matter. In a first 
approximation, hadrons should experience mass shifts, but considerable broadenings 
or strong modifications (e.g. mixing) of the spectral strengths are conceivable, too. 
A verification and understanding of such in-medium modifications of hadrons would 
allow to contribute to the "mystery of particle masses". 

This topic was essentially triggered by the Brown- Rho scaling hypothesis JTJ accord- 
ing to which the masses of the light vector mesons (V = p, lu, 0) behave as rriy oc (qq), 
where (qq) is the chiral condensate. In such a way fundamental and non-perturbative 
QCD quantities can be accessible. This idea influenced strongly the HADES project 
which aims at measuring the density and temperature dependence of (qq) via vector 
meson " spectral lines" . A strong impact on the field had QCD sum rules a la Shifman, 
Vainshtein, Zakharov [[J which, indeed, correlate condensates (and moments of the 
parton distributions in hadrons) and the pole masses of hadrons. 

The purpose of the present contribution is two-fold. First, we revisit the QCD sum 
rule approach to estimate the change of the vector meson pole mass parameter at finite 
baryon density n and temperature T. Second, we present dynamical calculations, bas- 
ing on hydrodynamics (fire ball model) and transport calculations (BUU) to elucidate 
the prospects to verify in the starting experiments at HADES the predicted in-medium 
modifications of vector mesons. 
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2 QCD sum rules revisited 



The basic object is the retarded current-current correlation function 

n*( g; ^,r) = i / d^e^e(x )([j^x) : j u (o)]) m , T , (i) 

with conserved vector currents = ^(ify^u =F d'j^d), S'j^s carrying the quan- 

tum numbers of the p, u and meson, respectively, but expressed by quark fields. 
(• ■ ')n N > T means grand canonical averaging with respect to chemical potential pn and 
temperature T. Here we focus on vector mesons at rest and constrain ourselves on the 
longitudinal part Ilf of the correlator IT?,,. 

The sum rule approach rests on a comparison of two different evaluations of Ilf . 
One can use (i) the analyticity properties to note a dispersion relation 



U^;, N ,T) = - ds '■r-" X :' ) (2) 



() , /M _ 1 r°° ds lmn*(s;n N ,T) 
it Jo s — (q + ie) 2 

where ■ • • indicate subtractions, and (ii) the Operator Product Expansion at large 
Q 2 = -ql > 

nf(Q 2 ) = -ainQ 2 + E^r, (3) 

n=l ^ 

where Cn are quantities containing Wilson coefficients and condensates and moments 
of parton distributions. Defining the spectral density phad(s; Pn, T) = ^Imllf (s; pn, T) 
and performing a Borel transformation one arrives at the QCD sum rule (QSR) 

ds p had (s) e-^ 2 = M 2 [C + £ ^h^ j > ( 4 ) 

where M is the Borel mass. Assuming that the r.h.s. is calculable from first principles 
one can use the QSR (i) either as consistency check of hadronic models for Phad(s) (cf. 
0), or (ii) as one equation for a physically relevant parameter in a simple parameter- 
ization of phad(s), or (iii) as constraint for a few parameters entering models of phad(s) 
(e.g. to find a width-mass relation ||, or to interrelate various parameters with the 
chiral gap |§]). 

We follow here the possibility (ii) as exercised by many previous authors like in 
0, ||. Inspired by the experimental results on the cross sections of the reactions 
e + e~ — > nn, K + K~ , which is intimately related to phad and which exhibits one low- 
lying distinctive resonance (either p or u or 0) followed by a flat continuum, one 
parameterizes 

PhadO; T > n) =F V 5(s- m 2 v ) + C Q(s - s%) + (p^t + Pscm)S(s), (5) 

where my is the wanted pole mass parameter, Sq stands for the continuum thresh- 
old, and p^cat^ denote the forward scattering amplitudes or Landau damping terms, 

„P,N _ /AQ^2\-l fOO 2 Al ,2 ft A 1\/T2 7, ,2 TO I A 1\/T2 /, ,21 



which read for the p meson p p s £ tt = (487H)" 1 du 2 h F Jl - AM 2 n /uj 2 [2 + AM 2 n /uj 



p££ tt = (247T 2 )- 1 IZl d^ 2 n B y/l - Am 2 Ju 2 [2 + Am 2 Juj 2 ] with n F = [e^~ 2 ^V 2T + 1]- 1 
h B = [e w / 2T — while for the meson, p^ltt — due to symmetry, and p^ t = 9Pseat 
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due to the different isospin structure of the p and u> mesons H; for the meson the 



Landau damping for the pion and nucleon gas is negligible ||10 



Truncating the r.h.s of eq. (f|) the sum rule might be cast in the form 



Co 1 - [1 + ^}e~ s °/ M - C 2 M 1 - C 3 M- 6 

m 2 = M 2 - - 

V C (l- e-oV^ 2 ) + C X M~* + C 2 M~* + |C 3 M-6 - + p^ tt )M~ 2 ' 

(6) 

where the coefficients up to mass dimension 6 have been calculated by many authors 
in the low-temperature or low-density approximation. In the low-T and low-n/y ap- 
proximation our result is 

C 2 = <?2 + 92 + 02, C 3 = q 4 + a 4 , (8) 

(12 



m q (uu} + 2M N a N YI» + ^m^Hj^ (9) 

g2 = I(V)o--M JV M fl Jf--^ (10) 
yi 24 7f 27 18 
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g4 = ira s Ko{uu) 



l _ k n 4M N a N yjN _ 36 

«« m,(-(ra) ) 1 7/2 1 



a 2 = ^^^+3^/^ + 4^^^ + ^^ ( 12 ) 

a 4 — — 2 JW 7Y /1 4 i AT 1 1 ^ N 4 > N 2 3 4 > JV 3 

^ i? " AZ d I* - 4A«+ d ij, (13) 

where ^ = 1, F = 1 for p (u) mesons, and elsewhere = 0. The sequence 
of replacements m q -> m s , (mm) -> (ss) K -> ym s /m q , and A^ >7r) -> A* i(iVi7r) 
holds for the <fi meson. The above integrals are I£ = J d 3 k[(2ir) 3 Ek\~ 1 k 2n ~ 2 uf, 1^ = 
J d?p[(2nYE P Y l p 2n - 2 n B , where n B = [e Ep/T - 1] _1 and n F = [e^-^/T + are 
thermal Boson and Fermion distributions, and /i/y is the chemical potential related to 
the nucleon density n N . 

Specific for our approach is the treatment of the four-quark condensate for which 
we extend the widely used ground state saturation assumption for the in-medium four- 
quark condensate at T = in the following way 



16 

((qi^ a q) 2 ), N = -{(qi,\ a qf), N = -K{n N )(qq)l N , (14) 



where the density dependent factor n{n^) is introduced to control a deviation from the 
exact ground state saturation. As pointed out in || n(n N ) reflects the contribution 
from the scalar low-energy excitations of the ground state and seems to be weekly 
dependent on n/y, so that one can use k(tin) = const as first approximation. In 
this case the density dependence of the four-quark condensate appears only via the 
density dependence of the quark condensate squared. In linear density approximation 
it is given by {qq) 2 m = (qq)o + (qq)o(N\qq\N)jj^. We go beyond this approximation 
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Figure 1: The in-medium changes of condensates as a function of temperature and den- 
sity (left: chiral condensate, top-right: strange chiral condensate, bottom-right: gluon 
condensate). Note that the validity range of the present calculations is constrained 
to small densities and temperatures, < 3n and T < 100 MeV are relevant for 
HADES. 



and include here a possible linear density dependence of n(n N ), keeping in mind that 
still «(njv) = k at = in accordance with the vacuum case. In leading order 
our extended ansatz results in an additional density dependence parameterized by a 
constant factor k,^. Our parameterization of the four-quark condensate at T = has 
then the form 

The limit = kq = 2.36 is used in ||, while the parameterization = 1.4 and 
k = 3.3 with (qq)o = (— 230MeV) 3 corresponds to the choice in f7j. Below we vary the 
poorly known parameter to estimate the contributions of the four-quark condensates 
to the QSR. The needed ansatz for the nucleon matrix element of the scalar four-quark 
condensate can be extracted from the direct comparison of our parameterization in 



eq. (|15|) and the general expression for the condensates via the matrix elements (the 
latter ones to be taken at T = 0) as (N\(qj^ 5 X a q) 2 \N) = ~(qq) (N\qq\N) k n . Since 
this matrix element does not depend on particle momenta and temperature it can be 
employed also for evaluating the four-quark condensates in the general case with T/0 
and n jv 7^ 0. 

The remaining parameters are specified by a s = 0.38, m q = 0.0055 GeV, m s = 
0.130 GeV, f n = 0.093 GeV, M N = 0.938_GeV, M% = 0.770 GeV, a N = 0.045 
GeV, y = 0.22, m n = 0.138 GeV, (uu) = (dd) = (-0.245 GeV) 3 , (ss)o = 0.8(uu) , 
(^G 2 ) = (0.33 GeV) 4 , A% + N d = 1.02, Af£ = 0.12, A u 2 % d = 0.97, A\% d = 0.255, 
A\ N = 0.1, A\ )N = 0.004, A| >7r = 0.08, A^J= 0.008. 

The density and temperature dependence of the chiral condensates (uu) and (ss) 
and of the gluon condensate (G 2 ) are exhibited in fig. 1. Note the linear density 
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Figure 2: The pole mass parameters m PjLJ ^ as a function of the density for T = 20 
MeV (solid curves) and 140 MeV (dashed curves) for kn/k = 0.42. 
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Figure 3: The pole mass parameter rn^j as a 
function of the density for = • • • 4 from 
top to bottom (T = 20 MeV, k = 3). 
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dependence and the nearly independence of the temperature. 

The parameter in eq. (|6|) is fixed by requiring maximum flatness of my(M) 
within a suitably chosen Borel window M min • • ■ M max determined by the " 10% + 50%" 



rule | j~lf . If we follow || and choose Kn/kq = 0.42 and use kq to adjust the vacuum 
values of my to the physical values, we get the results displayed in fig. 2. Due to some 
cancellations of the coefficients entering the p meson sum rule, the dropping mass 
parameter is determined by the four-quark condensate. For the u meson the Landau 
damping is dominating M. The dropping meson mass parameter reflects directly the 
strange chiral condensate; a smaller value of y diminishes strongly the mass shift. 

It turns out that the actual p and u mass shifts depend on the poorly known pa- 
rameter Kjv which parameterizes the density dependence of the four-quark condensate. 
In contrast to the above ratio kn/^o = 0.42 claimed in [||], the groundstate saturation 
hypothesis requires kn/^o — 1- This gives a stronger dropping m p . For the uo meson 
the rising m u , advocated in ||, is changed into a drop, as seen in fig. 3. Thus, even 
within the present simple pole approximation, the qualitative behavior of the u meson 
remains unsettled. 



3 Estimates of dilepton spectra in heavy-ion colli- 
sions 

Up to now we focused on the behavior of my within the QSR approach in zero-width 
approximation. The width-mass relation derived in is consistent with dropping 
p mass as long as the width is constrained to V < 600 MeV. However, a unique 
prediction of the in-medium modification of the width and the mass separately seems 
not be possible within the QSR approach. Rather, one has to resort to hadronic models 
which predict a considerable in-medium broadening. To elucidate the consequences of 
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Figure 4: Dilepton emission rate at T = 70 MeV as a function of invariant mass. Left 
panel: integrated over all M±, right panel: for M± = 1.3 GeV (dashed curve: v r = 0, 
solid curve: v r = 0.3). A vector meson broadening factor of 3 is assumed. 



these in-medium changes together we use the following parameterization of the emission 
rate of hadron matter [l^ 



dN 2d v f u-Q\l (m v T t - t )(MT v ^ e+e ^ 

d A x d 4 Q V (2vr) 3 6XP I T J vr (M 2 - m 2 v ) 2 + {m v Y^) 2 ' 1 ' 

where Q is the four-momentum vector of the lepton pair with invariant mass M, u 
denotes the four-velocity of matter, and dy stands for the degeneracy factor of the 
vector meson V. The corresponding invariant mass spectrum is displayed in fig. 4 for 
T = 70 MeV, where we assume that all vector mesons experience the overall broadening 
Ty 1 = bTy 1,0 with 6 = 3. As seen in fig. 4, the p meson is smeared out to such a degree 
that no distinct peak is visible. The u and can be identified as they are still sufficiently 
narrow. It is the double differential spectrum dN/d A x dM dMj_ which allows access to 
the p peak. As seen in the right panel of fig. 4, for fixed M± = 1.3 GeV there is a 
pronounced p peak with the to sitting on top. As known for a long time, such double 
differential spectra are sensitive also to the flow of matter. This is highlighted by a 
comparison of two flow velocities, v = (lower curve in left panel of fig. 4) and 0.3 
(upper curve). The full invariant mass spectra are insensitive to flow effects. 



3.1 Hydrodynamics 

In |12|] we have studied, within a scenario of an expanding fire ball, the resulting dilepton 
spectra. Instead of solving the hydrodynamical flow equations in full detail we have 
used a schematic expansion dynamics governed by R(t) = Ri + v r t with v r = 0.3, 
entropy conservation, s/n = const, and baryon conservation (N = 330). The final 
(freeze-out) parameters are from an analysis of hadron multiplicities at SIS energies 
[|T3|] nj = 0.3n , Tj = 50 MeV. Assuming an initial density rii = 3n , as suggested by 
BUU calculations, one gets Tj = 90 MeV for the initial temperature on the isentropic 
line. Keeping the above assumed in-medium broadening by a factor 6 = 3 and using 
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Figure 5: Space-time integrated dilepton spectra from the expanding fire ball described 
in text. Left panel: invariant mass spectrum integrated over all M±, right panel: 
M± = 0.9 GeV. A vector meson broadening factor of 3 is assumed and a p meson mass 
shift ofAm p = 300 MeV. 



the p meson mass shift of Am p = 300 MeV (a maximum value with respect to the 
above QSR studies) one gets the dilepton spectra exhibited in fig. 5. According to the 
above QSR results for kn ~ 2 we employ here no mass shift of the u meson. While in 
the invariant mass spectrum again no p peak is visible, the double-differential spectrum 
with Mj_ = 0.9 GeV shows clearly the shifted peak. 

In [12 also such scenarios are studied in detail, where the density depend broadening 
and mass shifts are incorporated. These studies support the need to measure double- 
differential cross sections to have access to the p peak. 



3.2 BUU calculations 

The above hydrodynamical calculations are based on local thermal and chemical equi- 
librium. These stringent assumptions are not longer needed in transport model cal- 
culations. We employ here a BUU code to calculate the dilepton spectra for central 
collisions Au(l AGeV) + Au. The results are displayed in fig. 6 for b = 3 and either 
no mass shift (left column) or a p mass shift of 150 MeV (right column). We have 
included the dileptons from 7r + 7r~ annihilations via the intermediate p meson, decays 
of p mesons from baryon resonance decays, pn bremsstrahlung, A and rj Dalitz de- 
cays and uj decays. In contrast to the above hydrodynamical calculations the p and u 
mesons are not allowed to undergo equilibration in the BUU code, but rather are prop- 
agated in mean field approximation along their original trajectories. The present BUU 
calculations confirm the results of the hydrodynamical model: Only when selecting a 
suitable M_l bin one gets access to the p peak. However, even in case of no mass shift, 
due to phase space weighting the p peak seems shifted (see left column in fig. 6). The 
role of the phase space weight is highlighted in the right column of fig. 6, where the p 
peak appears down-shifted by somewhat more than 150 MeV. 

In addition, the BUU calculations show that the Dalitz and bremsstrahlung yields 
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Figure 6: Dilepton spectra as a function of invariant mass for several M± bins. Left 
column: no mass shift, right column: density independent mass shift Am p = 150 MeV. 
An overall vector meson broadening by a factor 3 is assumed, u and contributions 
are not shown. do scal /dM = irR 2 (A projetile * At^)' 1 J AM± dM L dN/(dMdM ± ). 
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Figure 7: Dilepton spectra from 
uj decays as a function of invari- 
ant mass for several Q± bins. The 
density dependent mass and width 
are parameterized as = — 
8m w (n/n ) andT w = T^-5T w (n/n ) 
with Sm w = 70 MeV and ST^ = 
50 MeV. The sharp peak structure 
depicts the vacuum decay contribu- 
tion. 



do not disturb the wanted p signal. 

Various estimates demonstrate that the uj meson, even for some broadening, remains 
visible as distinctive sharp peak. Therefore, also small shifts should be visible. Here, 
however, a subtle interplay of uj decays in medium and after freeze-out occur. To 
elucidate the chances to verify possible in- medium uj modifications we have assumed 

= m" — 8mu(n/no) and = T° — ST^in/rio) with Sm^ = 70 MeV and ST^ = 50 
MeV and obtain the results exhibited in fig. 7. One can see clearly a double-peak 
structure from vacuum decays and in-medium decays. As expected, selecting low-Q^ 
dileptons the vacuum contribution is strongly reduced since the parent uj is nearly at 
rest. The smearing of the uj due to the density dependence can mimic a peak structure 
at the vacuum p position. 



4 Summary 

In summary we have revisited the QCD sum rule approach to the in-medium behavior 
of the light vector mesons in zero-width approximation. Due to cancellations the p 
meson mass shift depends on the four-quark condensate, but is always decreasing with 
increasing density. The large Landau damping and the uncertain four-quark condensate 
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do not allow a definite statement on the lo meson mass; for the ground state saturation 
the mass drops. The actual mass shift of the <fi meson depends sensitively on the 
uncertain strangeness fraction in the nucleon; the mass is directly related to the 
strange chiral condensate. 

Our hydrodynamical fire ball calculations show that a double-differential cross sec- 
tion needs to be measured to get access to the p peak in case of strong in-medium 
broadening. This is supported by transport calculations relying on a BUU code. To 
amplify a possible in-medium u shift the selection of low-momentum pairs is very 
helpful. 
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